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After reviewing the main characteristics of the spacetime of accelerating universes driven by a 
quintessence scalar field with constant equation of state to, we investigate in this paper the classical 
stability of such spaces to cosmological perturbations, particularizing in the case of a closed geometry 
and equation of state ui = —2/3. We conclude that this space is classically stable and conjecture 
that accelerating universes driven by quintessential fields have "no-hair". 
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I. INTRODUCTION 

Knowing the fate of the universe has always possibly 
been the main goal of cosmology. The speculation on 
what the future of the universe might be has now become 
even more complicated. Following recent observations in 
supernovas la [1] a real revolution has stirred cosmology, 
leading to a plethora of new views, concepts and finally 
to an emerging entirely new cosmological scenario which 
some like to call the New standard Cosmology [2]. Criti- 
cally pivotal of such a scenario is the already overwhelm- 
ingly accepted view that the expansion of the present 
universe is accelerating [3]. Whether this acceleration 
would continue forever or it will stop to recover deceler- 
ating expansion is still not known. The key ingredient 
that provides the repulsive force required by an acceler- 
ating universe is cosmic dark energy which, even though 
is not directly detectable, should make nearly the sev- 
enty percent of the total energy content in the universe. 
So far, two main candidates have been considered for 
dark energy: a positive cosmological constant [4] and the 
so-called quintessence field, a slowly-varying scalar field 
with negative pressure which only recently (in cosmic 
time) has started to dominate over all forms of matter 
[5,6]. Several compelling arguments have been advanced 
[5-7] in favor of cosmic quintessence with respect to the 
cosmological constant. 

It has recently been argued however [8-10] that if accel- 
erating expansion continues eternally, then a future event 
horizon will inexorably form and this would represent the 
demise for any mathematically consistent formulation of 
quantum gravity and string theory, both in the case that 
the universe be now dominated by a positive cosmological 
constant or by a cosmic quintessential field with constant 
[5] or time-dependent [6] equation of state. This possibly 
is one of the greatest challenges ever posed to theoretical 
physics. The formation of a future event horizon, which 
would make it impossible to construct any consistent S- 
matrix for string theory, is an unavoidable consequence 
from the causal future development of any initial sur- 
face in the accelerating universe if quantum coherence 
is to be preserved. Nevertheless, there exist solutions 
to the static Einstein equations for an eternally acceler- 
ating universe endowed with a quintessence scalar field 



with constant equation of state which do not show any 
event horizon [11]. This is made possible because such 
solutions possess a Kleinian signature that can lead to 
allowance of world lines traveling backward in time. Rel- 
ative to the case of a positive cosmological constant, this 
may become one of the strongest arguments in favor of 
quintessence. 

Along the evolution of the universe from its quantum 
era, tracking models can actually be represented [12] 
as a succession of uj = p/p= Const, domains, ranging 
from u> = +1 to uj = —1, separated by abrupt jumps 
and followed by the present accelerating expansion with 
— 1/3 > uj > —1, which evolves from an attractor solu- 
tion [6,12,13] and may or may not be characterized by a 
constant uj. However, it seems to be a good enough ap- 
proximation to also represent the present period of dark- 
energy dominated evolution by a quintessence model with 
w=Const. Therefore, investigating the spacetime of eter- 
nally accelerating universes whose expansion is driven 
by a quintessence scalar field with constant equation of 
state and does not show future event horizon appears to 
be of interest. This will be done in the present paper, 
specializing in the case uj = —2/3. We shall also deal 
with the issue of the stability of this eternally acceler- 
ating universe by considering the cosmological Lifshitz- 
Khalatnikov perturbations [14] on it. We check that the 
considered spacetime is stable and advance the conjec- 
ture that an eternally accelerating closed universe has 
no-hair. 

The paper can be outlined as follows. In Sec. II 
we review the solutions to the Einstein equations cor- 
responding to a quintessence scalar field with constant 
uj, minimally coupled to Hilbert-Einstcin gravity, both 
in the static and cosmological cases, for the whole range 
of state equations that covers all possible universes with 
accelerating expansion. We particularize in the case 
uj = —2/3 whose FRW metric is considered in some detail 
for the different geometries of the universe. The Lifshitz- 
Khalatnikov cosmic perturbations of this spacetime are 
studied in detail in Sec. Ill and it is checked that the 
spacetime is stable to them. We also advance the con- 
jecture that an eternally accelerating closed universe has 
no-hair, discussing it by comparing with the purely de 
Sitter space. Finally we conclude and add some further 
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remarks in Sec. IV. 



II. THE SPACETIME OF AN ACCELERATING 
UNIVERSE 

The spacetime of the accelerating Friedmann- 
Robertson- Walker (FRW) universe endowed with a 
quintessence scalar field with constant state equation pa- 
rameter — 1/3 > cj > —1 corresponds to maximally sym- 
metric spaces with negative spacetime curvature which 
are solutions of the Einstein equations. Since the conser- 
vation law for a generic quintessence scalar field should 
be taken to be ap = a(t)~ 3 ( 1+0J ^ (where a is an arbi- 
trary integration constant, p is the energy density of the 
quintessence field and a(t) is the time-dependent scale 
factor), the Friedmann equations for the quintessential 
accelerating universe are [15] 
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3a 



3a 



-3(l+w) 



= 



± 



R 2 a 2 



(2.1) 



(2.2) 



in which H is the Hubble constant, R~ 2 is the spatial 
curvature constant, and the overhead dot denotes dif- 
ferentiation with respect to time t. Note that if we set 
uj = — 1 the first term of the right-hand-side in Eq. (2.2) 
becomes a constant, and we obtain then a solution for 
the scale factor that describes just de Sitter space. Here 
we shall restrict ourselves to solve Eqs. (2.1) and (2.2) in 
two particular interesting cases. When we set u) = —1/3, 
i.e. at the onset of the accelerating regime, we obtain the 
solution [16] 
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t + K , 



(2.3) 



where K is an integration constant. Whereas in the 
spatially closed case with i?~ 2 = 8ttG/ (3a) the scale fac- 
tor (2.3) reduces to a simple constant that describes an 
Einstein static universe, in the spatially flat and open 
cases, or when i?~ 2 < 8irG/ (3a) for closed geometry, the 
universe will expand in just the uniform way. 

For a constant equation of state with lu = —2/3, i.e. at 
the typical most interesting situation in which the uni- 
verse expands in an accelerating fashion quite adjustable 
to what has been observed in recent supernova experi- 
ments [1,3], the scale factor solution to the Friedmann 
equations (2.1) and (2.2) can generally be written 



, , 2irGt 2 rjr 3a (K 2 + iT 2 ) 
a{t) = — + Kt + 
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(2.4) 



If, without any loss of generality, we set the integration 
constant K = (in what follows it will be seen that when 
we re-express the scale factor in terms of the compactified 
time 77, the resulting solution does not explicitly depend 



on K), then this solution described the FRW spacetime 
of half a Lorcntzian wormhole [17], from a throat at t = 0, 
where a = clq = 3a/ (8irGR 2 ), to the asymptotic region 
at t = +oo. We furthermore notice that if for K = 
we allowed time t to take also on negative values from 
t = down to t = — oo, then a complete wormhole would 
be obtained. Let us next assume for a moment that the 
two asymptotic regions at t = ±oo are identified to each 
other so that [18] any world lines for test particles or light 
signals approaching the infinity surface at t = +oo from 
t = would find themselves emerging from the infinity 
surface at t = -co, toward t = 0, and any of such lines 
approaching the infinity surface at t = -co also from 
t = would emerge from the infinity surface at t = +co, 
again toward t = 0. None of these world lines had then 
reached any of the infinities. Thus, if the universe can be 
described as a compactified complete wormhole the way 
we have just described, then light signals traveling back- 
ward in time existed and could connect the whole space- 
time in such a way that no event horizon would form up 
in the future [11]. This is by no way implying the exis- 
tence in the accelerating universe of any nonchronal re- 
gions containing closed timelike curves because for these 
curves to occur in our spacetime it would be necessary 
that the two asymptotic regions at t = ±oo were also 
set into motion relative to one another [18], which is not 
possible for the universe. 

In what follows, we shall interpret the throat of the 
wormhole at t — as the latest hypersurface of the uni- 
verse immediately before the onset of the accelerating 
regime. This will make the value of the constant 3a/87rG 
very large. The three possible geometries associated with 
the FRW metric that correspond to solution (2.4) with 
K 7^ lead to the following expressions for the scale 
factor in terms of the conformal time r) = J dt/a(t). 
(i) Spatially closed R~ 2 > 



a(v) = 



3a 



8irGR 2 cos 2 (ry/(2i?)) : 



with 



r] = 2R arctan 



R(K + 



4wGt \ 
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(ii) Spatially flat R- 2 = 



with 
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rj = - 



K + 
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(iii) Spatially open R 2 < 
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Note that in all three cases the scale factor does not ex- 
plicitly depend on the integration constant K. 

Although particular values of the parameter lu such as 
lu = —1/3 and u> = —2/3 would strictly make the usual 
sense only for homogeneous and isotropic FRW space- 
times, if we keep an equation of state p — Lup also in 
the case of spacetimes with static, spherically symmetric 
coordinates, one can also obtain the static metrics which 
correspond to the above particular cases for a given, fixed 
relation between the energy density and the metric com- 
ponents, much in the same way as the static metric for 
de Sitter space can be derived from the Einstein equa- 
tions for static, spherically symmetric coordinates and 
an equation of state p = up, whenever we set to = — 1. 
That static de Sitter metric can directly be related with 
a corresponding FRW metric for de Sitter space obtained 
from Eqs. (2.1) and (2.2) also for 10 = — 1 by means of 
embeddings in a common five-dimensional hyperboloid 
(see Refs. [20] and [21]). We next consider the generic 
metric that describes the spacetime in static, spherically 
symmetric coordinates t, r, 8, cf) for an equation of state 
p = ujp and positive spatial curvature, corresponding to 
an ansatz 

ds 2 = -A(r)dt 2 + B(r)dr 2 + r 2 dfll, 

where dfl 2 is the metric on the unit two-sphere, and A(r) 
and B(r) are metrical coefficients depending only on the 
radial coordinate r. Following Rcf. [11], we take as the 
components of the energy-momentum tensor T,!' = Tg = 
T$ = p and T ° = -p, that is T r = T e e = = -wT§, 
and hence we have the Einstein equations 



1/1 A' 
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in which the prime denotes derivative with respect to 
the radial coordinate r. On the other hand, the compo- 
nents of the energy-momentum tensor must satisfy the 
gravitational equation [11] T* k = 0. From this equa- 
tion wc can now obtain a relation between the met- 
ric component A(r) and the energy density p, A(r)' = 
—2u;A(r)p'/[p(l + lu)], which can be immediately inte- 
grated to give ap = A{r)~^ + ^^ 2oj \ with a again an 
arbitrary integration constant. Note that for to = —1, p 
becomes a simple constant, so that the Einstein equations 
are straightforwardly solved to produce the well-known 



static de Sitter metric. In the case co = —1/3 we can 
obtain the solution [11,16] 
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where K\, Ki and K% are integration constant which 
are arbitrary unless by the condition KiK 3 = —1. This 
metric shows an event horizon at 

_ (3aK 3 K 2 8irG V 1 
r ~ Th ~ { 8ttGK 2 ~ 3a~K^ ) ' 

which marks the transition from a Kleinian-signature 
metric for r < ru to an Euclidean (Riemannian) met- 
ric for r > rh. For the entire range of lu- values in the 
accelerating interval —1/3 > lu > —1, we have the solu- 
tion [11] 



B 



u> 2 + 6lu + 1 

(1*LU) 2 ' 



where 
K{lu,o) 



f27rG(l + Lu) 2 B\ 2uj/(1+uj) 

— , -Kw<-l/3 

\ ua J 



Note that in the considered interval B reduces to a sim- 
ple dimensionless constant which is negative definite. Be- 
sides being singular at the origin of radial coordinate r, 
the two most interesting properties of this solution are: 
(i) it does not show any event horizon (so that, all world 
lines should necessarily be always connected to each other 
in both the static spacetime and, in spite of correspond- 
ing to an accelerating FRW spacetime, the cosmological 
spacetime described by metric (2.4)), and (ii) it has a 
Klcinian definite signature ( h +). The latter prop- 
erty is consistent with property (i) and with a Lorentzian 
wormholc interpretation of the FRW metric with scale 
factor (2.4) for the following reason. Let us first in- 
troduce the coordinate change r = £t&n 2 (ip/2), where 
I = K(lu, a) _ ( 1+w )/( 4w \ with the new half-periodic coor- 
dinate V> running from = 0(r = 0) to = ir(r = oo), 
and then consider a world line for matter defined by 
8, </>=const., t = —flip, in which /3 is a real constant. The 
static metric reduces then to the line element 



ds 2 = - 



(3 2 tan : 



8w/(l+w) 



(V>/2) - Bt 



tan 2 (V>/2) 
cos 4 (W2) 



d^ 2 



Since B is a negative definite constant along the entire 
interval —1 < lu < —1/3, and (3 is real, this line element 
is always timelike along that interval. Thus, an observer 
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moving on the world line will always have an increas- 
ingly negative time coordinate, i.e. though that observer 
cannot even reach the point ip = ir and, therefore, can 
never follow a closed timelike curve, it will always travel 
backward in time. 

In this paper, we shall concentrate on the particular 
value to = —2/3 which corresponds to a FRW metric that 
predicts an accelerating universe which conforms well to 
supernova results. In this particular case, the above so- 
lution reduces to 

which retains all the properties which we have discussed 
for the generic case. 

Most recent results on the large scale curvature of the 
universe provided by BOOMERanG [19] and other CMB 
anisotropy experiments [20] indicate that the universe is 
fiat with a 95 percent confidence interval. This leaves a 
5 percent uncertainty which, from the qualitative stand- 
point, tells us that the geometry of the universe can still 
be open or closed. On the other hand, along the devel- 
opment of physics we are used to finally discover that 
the apparently most probable simplest situations (circu- 
lar orbits for planets and atomic electrons, spherical sym- 
metry, etc) were often not the real case, but just a good 
approximation. If flatness is taken to represent the sim- 
plest possible geometry of the universe, we shall adhere 
to follow this empirical tendency by choosing for our ac- 
celerating universe with quintessential equation of state 
lu = —2/3 a closed geometry, hoping that, like for pre- 
vious historical examples, this finally uncovers a richer 
structure. Moreover, one should expect that the conclu- 
sions on the stability of the universe against cosmological 
perturbations we are going to obtain in the next section 
for closed geometry would be also shared by flat geom- 
etry In any event, a detailed consideration of the cases 
for flat and open geometries are left for future research. 
It is only for the sake of completeness that we shall fi- 
nally include the metrical solution for any constant to in 
the accelerating interval for the flat geometry. In FRW 
coordinates it reads 

ds 2 = -dt 2 + I i '- I (dr 2 + r 2 dn 2 2 ) . 

(2.14) 

This line element is valid for all —1/3 > to > — 1 and is 
associated with a corresponding static metric for any ac- 
celerating value of co [11] of which the metric coefficients 
in Eq. (2.13) become just the case for to = —2/3. There- 
fore, metric (2.14) and its static counterpart should pos- 
sess exactly the same properties as what were discussed 
before for metrics (2.4) and (2.13). 



III. COSMOLOGICAL PERTURBATIONS 

In this section we shall first briefly review those math- 
ematical aspects of the perturbative technique developed 
first by Lifshitz [14] which will be later used to investi- 
gate the stability of our closed accelerating universe. For 
a closed FRW spacetime, perturbations of the four-metric 
are introduced by g a b — gab + h a b, with the perturbations 
satisfying the gauge hoo = = 0, a,b, ... = 0,1,..., 3 
and a, (3, ... = 1,...,3. It is work remarking that, al- 
though we are working in an isotropic framework with a 
reference system which is always synchronous, after this 
choice of gauge, it is no longer a comoving system be- 
cause the perturbations on the spatial components of the 
four- velocity can be generally nonzero [14,22,23]. Pertur- 
bation of the Ricci tensor and, hence the Einstein equa- 
tion, energy density and velocity can then be derived (for 
details see Refs. [21,22]). By taking as the most general 
metric perturbation 

h a p = Hv) P al3 + V(v)Qa0 + Cr(v)Sa/3 + v(rj)H a p, 

(3.1) 

where the coefficients A, fi, a and v depend only on the rj 
time parameter, and P a /3, Qa/3, S a p and H a p are tensor 
harmonics which arc defined by [22] 

Q a p = \lapQ, V a V a Q(") = (l-n 2 )Q(") (3.2) 

Pa/3 — + Qa/3, S a f} — V a S ' f) + V pSa, (3.3) 

V Q V a ^ n) = (2 - n 2 ) S { b n \ V a 5i n ' = (3.4) 

V„V»^;) = {3-n 2 )H%\ V a ^=0, H^ a = 0, 

(3.5) 

with Q^ n \ sl™' and , respectively, the scalar, vector 
and tensor harmonics [21,22], an n an integer order. 

Inserting Eq. (3.1) into the perturbed expressions for 
the Einstein equations and the energy density and ve- 
locity components, and using the above definitions, we 
finally obtain (a prime denotes differentiation with re- 
spect to rj) 

1n'\' 1 

A" + -e{l + 2)(X + n) = (3.6) 

ci o 

H" + {3C 2 + 2) ^ 

+ 1 (3C S 2 + 1) \t{i + 2) - 3] (A + M ) = 0, (3.7) 
a „ + 2a^ =Q (3g) 
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(3.9) 



and 



^^^{^ + 2)-3](A + M) + 3^0 
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pa 
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(3.10) 



{£{£ + 2)n' + [^ + 2) - 3] A'} , 



for scalar harmonics, and 
Sp 



= 0, Sv a = [£(£ + 2) - 3] a'S° 



(3.11) 



(3.12) 



for vector harmonics. 



A. scalar harmonics 



A = ^ + 2)(77 + sinr7), p, = -i 



+ 2) (i] + sin rj) + 3 sin r\. 

(3.19) 



These solutions will be conveniently subtracted from the 
general solutions to Eqs. (3.13) and (3.14). The latter 
solutions will be obtained with the help of the auxiliary 
functions £ and £ introduced by the substitutions 



X + p 



— 3 sin r] J 



(3.20) 



A'-// = 2^ + 2) (1 + cos 77) / £dr] -3cosry / £dr] + (cost). 

(3.21) 

When substitutions (3.20) and (3.21) are introduced in 
Eqs. (3.13) and (3.14), we obtain the new coupled differ- 
ential equations 



l--scc 2 (^/2) 



C = (3.22) 



The differential equations for the 77-dependent coeffi- 
cients A and /x, representing metrical scalar perturbations 
of a closed FRW geometry for the scale factor (2.5) can 
be written as 



A" + 2 tan(r?/2)A' - -£{£ + 2) (A + p) = 



-[£(£ + 2)-3] (A + M ) = 0, 



(3.13) 



(3.14) 



while the perturbations for energy density and velocity 
components become 



Sp COS 4 (77/2) 



5v c 



12 [1 - tan 2 (7?/2)] 



{[£(£ + 2) - 3] (A + it) + 3 tan W '} Q 

(3.15) 



{£(£ + 2)// + [£{£ + 2) - 3] A'} : 



(3.16) 

where we have consistently taken for the speed of sound 

C s = Sp/8p = uj = -2/3, (3.17) 

which corresponds to a quintessence scalar field with con- 
stant equation of state u = —2/3 [12]. 

The differential equations (3.13) and (3.14) still con- 
tain some residual gauge freedom for a complete speci- 
fication of the choice of coordinates [14,21,22]. Such an 
unphysical gauge corresponds in the present case to the 
particular solutions 



A = — p = Const 



(3.18) 



C + [tan(?]/2) - tan 77] C 

+ [2£(£ + 2) (1 + sec rj) - 3 (1 - tanr/tan(r//2))] 
= 0. (3.23) 

Straightforward manipulations on these equations lead 
finally to 



( = — sec r/ cot 2 {f}/2)y' 



cot(»?/2y 



2 -£{£ + 2)-\ + \t^{r i /2) 



(3.24) 

y = o, 

(3.25) 



where y = £ sin 2 (rj/2). 

Our task now is to solve Eqs. (3.24) and (3.25). Since 
obtaining exact solutions to these equations in closed 
form is very difficult, we shall derive approximate so- 
lutions in the extreme cases when 77 — > (i.e. at the 
beginning of the accelerating phase) and rj — > n (i.e. to- 
ward the asymptotic future of the eternally accelerating 
expansion). In the first stages of accelerating expansion 
with rj « 1, Eq. (3.25) can be approximated up to 
second order in rj as 



-rjy" + - (12 - rj 2 ) y' + 



0. 

(3.26) 



At the smallest rj the solution can in turn be approxi- 
mated in terms of Bcsscl function J v in the form 



y^rj 3/2 J: 



3/2 



(3.27) 
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As 77 — > 0, we then have y oc jf/ 2 . It follows that for the 
auxiliary functions £ oc 77, ( oc —4(1 — ?] 2 /3), and hence, 
since all residual gauge given by Eqs. (3.10) will vanish 
as 77 — > 0, one consistently concludes that the metric and 
dark energy (density and velocity components) perturba- 
tions all vanish as one approaches the onset of the accel- 
erating region. Had we chosen for the Bessel function any 
of the functions H v [23] , then all the above perturbations 
would be pure imaginary and divergent as 77 — > 0. 

Of greater interest to study the stability of our eter- 
nally accelerating universe is to consider the behaviour 
of perturbations in the asymptotic region r/ — > 7r (i.e. 
t — > 00). Thus, we next look at the solutions to Eqs. 
(3.24) and (3.25) as 77 — > w. Let us first introduce the 
change of time coordinate x = tan(7j/2) in Eq. (3.25) 
which then becomes 



-(! + ; 



!\2 <Py_ 
' dx 2 



x v 'ax 



-£(£ + 2) - 1 + 3x 2 



y = 0, (3.28) 



which for large x and even moderate £ admits an approx- 
imate solution again expressible in terms of the Bessel 
function J„, i.e. 



y/8£{£ + 2)- lcot(»7/2) 



V- v / cot(?7/2)J v ^3 /2 
As 77 — > 7r this solution reduces to: 

yoc[cot(77/2)] (1+v ^ )/2 



(3.29) 



or when expressed in terms of the auxiliary functions £ 
and £, 

£ ~ sin- 2 (77/2) [cot(7//2)] (1+y " )/2 (3.30) 



C ^ isccT? (l + >/l3) £, (3.31) 



where A(f) is a finite constant that depends only on £ and 
whose precise value is not of interest in this paper. The 
use of these expressions in Eqs. (3.20) and (3.21) leads 
finally to the solutions for the perturbations coefficients 



A = [B(£) - 12sin 2 (77/2)] [cot(r7/2)f+^ 2 

(3.32) 



„ = ~-^j= [B{£) + 12sin 2 (77/2)] [cot(?7/2)]( 5 +^/ 2 , 
o + v lo 

(3.33) 



where B{£) is a finite constant given by 



B{£) = 



16 



1-^13/8-^ + 2) , 
5 + v 13 L J 
and we have subtracted the unphysical gauge as- 
sociated with the particular solution given by Eqs. 
(3.19), i.e. A = 2C\£(£ + 2)(r) + sin??) + C 2 , [i = 
-2Ci [<?(£ + 2) (77 + sin 77) -3 sin 77] - C 2 , with d and 
C2 two arbitrary integration constants. Using then 
Eqs. (3.15) and (3.16) we finally derive the perturba- 
tion in energy density and velocity components for the 
quintessence field which are given by 



Sp _ 2A(£)cos 4 (r]/2) [cot(7j/2)] 



(5+vT3)/2 



3(3 + vT3)csc 2 (77/2) 

-4 [£(£ + 2) -3] - 12sin7;cot(77/2)+ 5 + [B(£) csc 2 (r]/2) + 12] 1 Q 



(3.34) 



Sv ° = m_cos 2 (r,/2) [cot(77/2)] (3+ ^/ 2 



[£(£ + 2) - 3] 



(3 + V13)C0S77 



+ 2) 



(5 + V13) 



(B{£) esc 2 (77/2) + 12) - 12sin77[cot(?7/2)] 



(2+^13)/2 



(5 + 4 Vl3) (B(£) csc 2 (77/2) + 12) + 12 sin 77 [cot(77/2)] (2+/T5)/2 



pa 



(3.35) 



All of the Eqs. (3.32) - (3.35) vanish as 77 — > w. Had we taken any of the two Bessel functions H v [23] instead 
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of J u for the solution of the differential equation (3.28), 
then we would have finally obtained the same expres- 
sions as (3.32) - (3.35), but with the sign for all \/l3 
changed. These would again vanish as r\ — > n, so the 
accelerating closed universe resulting from the presence 
of a quintessence scalar field with constant equation of 
state u = —2/3 appears to be stable to scalar Lifshitz- 
Khalatnikov perturbations. 



B. vector and tensor harmonics 

For the case under consideration, perturbations associ- 
ated with vector harmonics are described by coefficients 
that satisfy the differential equation 



cr" + 2tan(r?/2)cr' = 0, 



(3.36) 



with 



^ =0, 5v a = [l(i + 2)-'i]a'S a . (3.37) 
P 

The solution to Eq. (3.36) can be given in closed form 
and reads 

a = C + Ci ^r? + sinr? + | sin(2r7)^ , (3.38) 

where Co and C\ are arbitrary integration constants. 
Now, from Eqs. (3.37) it follows that a' = 0, so the con- 
stant C\ should be zero too, and hence a = Co- There- 
fore, as usual [14,21,22], vector perturbations correspond 
to unphysical pure gauge and are hence irrelevant also 
for the eternally accelerating universe we are considering 
in this paper. 

Let us next deal with the quite more interesting study 
of the gravitational waves associated with the perturba- 
tions generated by tensor harmonics. In our case, the 
coefficients for such perturbations are described by the 
differential equation 



v" + 2 tan(r?/2)z/ + £(£ + 2)u = 0. 



(3.39) 



Two cases can now be distinguished. If I = 0, then the so- 
lution for v(ti) is formally identical to that for <j(r]) given 
by Eq. (3.38), but its interpretation is different. It physi- 
cally represents the time evolution of gravitational waves. 
We have obtained that, even though the gravitational 
wave amplitude does not vanish as r\ — > w, neither it grow 
as t — > oo, at which limit v — C = C + 37rCi/4=Const. 
The neat effect of the whole evolution from rj = to 
rj = 7ron the zero-mode amplitude is an increase from 
Co to C*. If £ ^ 0, then by introducing the coordinate 
change z — sm(t]/2), Eq. (3.39) can written as 



(1 - z 2 )v" + 3zv' + m(m + A)v = 0, 



(3.40) 



in which m — 2£. The solution to this differential equa- 
tion can most easily be expressed in terms of ultrasphcr- 
ical (Gegenbauer) polynomials C'n ^ [23] of odd degree n 



and reads 

v( V ) = cos 5 / 2 (77/2)^ [sinfa/2)] , £ > 0. (3.41) 

We can then show that this solution becomes propor- 
tional to i~ 5 / 2 as one approaches the asymptotic limit 
r] = it. Thus, these gravitational modes (which all start 
with vanishing amplitude at rj = 0) become asymp- 
totically suppressed as one goes to t = oo. We have 
thereby excluded any unstable growing modes of the 
gravitational radiation in the accelerating regime driven 
by a quintessence field with constant equation of state 
u = -2/3. 



C. No-Hair conjecture 

Now, from the solution (3.38) for £ = we obtain at 
large t, 



-y (0) 
w=-2/3 



27a 
1287rGa 



Ci. 



(3.42) 



Thus, even though the value of v provided by Eq. (3.30) 
would be expected to be larger than the corresponding 
value for the zero-mode, = C 2 =Const., for de Sitter 
space [21], one may always choose the value of the con- 
stant C*,Ci and C2 such that ^f,°2_ 2 /3 became smaller 



than v*fo at sufficiently later times. For I ^ 0, from 
solution (3.41) at large t, we also obtain 



w=-2/3 



(21 + 4)! 
51(2^-1)! 



(2a) 



-5/4 



(3.43) 



The comparison of this with the corresponding de Sitter 
expression [21] ~ A\ + A 2 cxp(— 3Ht) (where H is 
the Hubble constant and A\ and A 2 are constants which 



-2/3 



< v. 



only depend on £) clearly implies that v t 
even at moderate values of time t. 

Moreover, from the above discussion it follows that 
any physical effects driven by the gravitational radia- 
tion modes (3.42) and (3.43) should in any event be very 
small, since their physical wavelengths respectively in- 
crease with a 1 / 6 and a 5 / 12 . Thus physical quantities in- 
volving at least two derivatives of the metric are then sup- 
pressed asymptotically by powers of the inverse of these 
wavelengths. These results appears to be implying a cos- 
mic "no-hair" theorem [24] for our accelerating closed 
universe endowed with a quintessence field with constant 
equation of state u> = —2/3. According, furthermore, to 
our discussion above, it would rather be a question on 
the relative values of the constants involved in the ex- 
pressions for coefficients v whether the no-hair de Sitter 
attractor or our no-hair attractor is the final solution for 
an accelerating universe. Of course, the model discussed 
in this paper is a simple one, so that one should extend 
this discussion to include other models with different spa- 
tial geometries and constant or "tracking" equations of 
state. 
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IV. CONCLUSIONS 

Motivated by the huge impact produced by supernova 
observations [1] , in this paper we have studied some new 
characteristics of the cosmological and static spacetimes 
of accelerating universes whose expansion is driven by 
a slowly-varying quintessence scalar field with constant 
equation of state —1/3 > u > — 1. Particularizing at 
the observationally most favored case u) — —2/3, it was 
shown that these spacetimes do not possess any future 
event horizon, a property which is not obviously shared 
by de Sitter space. The reason that justifies this result 
is that, even though closed timelike curves are not per- 
mitted to exist in these spacetimes, world lines traveling 
backward in time are allowed to occur on it, so connecting 
any otherwise causally disconnected regions, and hence 
preventing the formation of future event horizons. This 
result could have considerable interest for particle physics 
and quantum gravity because it manifestly avoids the 
serious, perhaps fatal difficulties for string theory (and 
actually any quantum field theory that depends on the 
presence of particles at infinity) posed by the existence 
of a future event horizon [8-10]. This would be just an- 
other reason in favor of preferring quintessence over a 
positive cosmological constant, since asymptotically de 
Sitter space inexorably leads to the formation of an event 
horizon. 

The stability of the closed space with accelerating ex- 
pansion induced by a quintessence field with an equa- 
tion of state u> = —2/3 to the cosmological Lifshitz- 
Khalatnikov perturbations [14] on the three-sphere has 
been also studied in detail. Although there exist more 
elaborated, covariant methods for dealing with cosmolog- 
ical perturbations [25], we have followed here the origi- 
nal Lifshitz-Khalatnikov treatment because of its greater 
adequacy to distinguish among the involved physical ef- 
fects. It was obtained that, at least for a closed geome- 
try in the case u = —2/3 the space is stable to both, the 



scalar and tensorial perturbations, and that the damping 
of small physical effects induced by the resulting gravi- 
tational waves allows one to conjecture that -much like 
it happens in asymptotic de Sitter space- eternally ac- 
celerating universes induced by quintessential fields have 
no-hair, so becoming final attractors along the evolution 
of the universe. It would be the value taken on by the 
constants that characterize the amplitude of the gravi- 
tational radiation what would finally decide whether the 
attractor of de Sitter universe or the attractors of lo > — 1 
quintessential accelerating universes are going finally to 
dominate and drive the future cosmological evolution. 
This is a matter which cannot be decided in the present 
paper, but that appears to be decisive to avoid the above- 
mentioned severe conflict between accelerating universe 
and string theory (or possibly any competing quantum 
theory, if there were any) . Since our present understand- 
ing of any of the theories involved at this conflict is still 
too rudimentary, it is still rather premature to say any- 
thing definite about it. 

It is expected that the results obtained in this paper for 
the closed geometry can be applied to flat geometry too. 
In particular, it appears that the no-hair conjecture be 
a key ingredient also for the spatially flat case. Finally, 
the rather intriguing implication that in any eternally 
accelerating universe driven by quintessence there would 
be world lines (followed by light signals and possibly some 
kind of matter) traveling backward in time need further 
consideration. After all, if such lines were allowed to 
exist, then the meaning of cosmological evolution itself 
should actually require a deep revision. 
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